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Abstract. By comparing the results from numerical microlensing simulations to the observed long-term variability of quasars, 
strong upper limits on the cosmological density of compact objects in the 10""* Afg-1 Mg range may in principle be imposed. 
Here, this method is generalized from the Einstein-de Sitter universe to the currently favored flu = 0.3, CIa = 0.7 cosmology 
and applied to the latest observational samples. We show that the use of high-redshift quasars from variability-selected samples 
has the potential to substantially improve current constraints on compact objects in this mass range. We also investigate to 
what extent the upper limits on such hypothetical dark matter populations are affected by assumptions concerning the size of 
the optical continuum-emitting region of quasars and the velocity dispersion of compact objects. We find that mainly due to 
uncertainties in the typical value of the source size, cosmologically significant populations of compact objects cannot safely be 
ruled out with this method at the present time. 
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1. Introduction 

Despite much effort to detect or rule out cosmologically sig- 
nificant populations of compact bodies in the stellar to plane- 
tary mass range, such objects still remain viable candidates for 
the dark matter of the universe. Although many of the com- 
pact dark matter candidates proposed are baryonic (e.g. red, 
white and brown dwarfs, neutron stars, stellar black holes, 
gas clouds) and therefore constrained by standard Big Bang 
nucleosynthesis and the determination of the primeval deu- 
terium abundance to contribute no more than Qb/iL) = 0-02 



(e.g. Buries et al. 2001) to the cosmological density, several 
candidates do circumvent these constraints, either by being 
non-baryonic or by making their baryonic content unavailable 
by the time of nucleosynthesis: e.g. primordial black holes 



(Hawking 1971 1, quark nuggets (Alam et al. 1999), mirror 



matter MACHOs (Mohap atra 1999) and aggregates of bosons 
or fermions (Membrado 1998 ). A large number of meth- 



ods do however exist to constrain the cos mologi cal densities 
of such populations (see Dalcanton et al. 
Sakellariadou 1999 for reviews). 



1994 and Carr & 



In this paper, we will be concerned with compact objects 
in the potentially very interesting mass range 10 Mq-1 Mq, 
where indirect detections of cosmologically signi ficant popu- 
lations have been suggested (e.g. Hawkins 1996 ). In this re- 
gion, the currently most powerful constraints on the cosmolog- 
ical density of such objects (regardless of type), ficompact, come 
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from Dalcanton et al. (1994) and Schneider (1993, hereafter 



S93), and are both based on theoretically predicted effects of 
quasar microlensing. 

By comparing the equivalent width distribution of quasar 
emission lines predicted by microlensing scenarios to that of 
observed samples, Dalcanton concludes Qcompact ^ 0.2 for 
lO"-' Mq-1 Mq compact objects in a critical universe. Taken at 
face value, the constraints from S93 are however even stronger 



at the lowest masses: Q 



compact 



< 0.1 for 10"' Mo-3 ■ 10"^ Mr.. 



The method of S93 is based on the argument that large pop- 
ulations of compact objects should statistically induce varia- 
tions in quasar light curves larger than those actually observed. 
The limits derived this way do however rely on the premise 
that the many parameters going into the microlensing simula- 
tions can be sufficiently well constrained by observations or 
reasonable assumptions, effectively making Qcompact the only 
free parameter The aim of this paper is to investigate whether 
this may actually be accomplished at the present time. Of par- 
ticular importance is the ill-determined typical size of the UV- 
optical continuum-emitting region of quasars, Rqso, which S93 
assumes to be Rqso = 10'^ m when deriving the constraints 
quoted above. 

Here, the method of S93 will be generalized from the 
Einstein-de Sitter (EdS) universe to the cuiTently favored Qm = 
0.3, Qa = 0.7 cosmology and applied to the latest observa- 
tional samples (Hawkins 200C| , hereafter H2000). The sensitiv- 
ity of the constraints to uncertainties in Rqso and the velocity 
dispersion of compact objects will be evaluated using recently 
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developed methods to better approximate the magnification in 
the case of large source microlensing. 

2. Numerical simulations of quasar microlensing 

2.1. Computational method 

The method used to derive the statistical properties of light 
curves of quasars microlensed by a cosmological distribution 
of compact objects is based on the machinery outlined in S93, 
extended to arbitrary Friedmann-Lemaitre cosmologies using 
the angular size distances of Kayser et al. ( 1997 ) and to the 
case of large-source microlensing using the magnification for- 
mula of Surpi et al. (2002). 

In this technique, the multipl icative magnification approx- 
imation (Ostriker & Vietri 1983 ) is assumed to adequately re- 
produce the statistical probability of variability. In this case, the 
magnification, yUtot, due to / microlenses is equal to the product 
of the individual ones: 



/^tot 



(1) 



Although this approximation cannot be expected to repro- 
duce detailed features of microlensing light curves, it has 
been proved useful (e.g. Pei 1993 1 for statistical investigations 
and significantly reduces the calculation time required by ray- 
tracing and similar methods. 

We follow conventions commonly found in gravitational 
lensing literature and use D^, Ds and Dds to denote the angu- 
lar size distances from observer to lens, observer to source and 
from lens to source, respectively. If ^ is the separation of the 
lens to the line that joins the source and observer, the dimen- 
sionless impact parameter y (in units of the Einstein radius of 
the lens) then becomes 



(2) 



where Mcompact is the mass of the compact object acting as a 
microlens. The dimensionless source radius R may similarly 
be expressed as 



R - cRq5o 



(3) 



For a circular source with uniform brightness, the magnifi- 
cation due to a single microlens is in the limit of small sources 



(R < 1.2) approximated by (Schneider 1993) 
2(R^+2) 



R(R^ + 4)3/2 



1 + 



(R^ + 2)/ 



R\R^ + 3) 



(R"-+3)^+3-4Cv/«)- 

«(4+R2)3/2 

2+/ 



fory < R 
for V > 



(4) 



In the limit of large sources (R > 1 .2) the following approxi- 
mation (Surpi et al. 2002 ) is instead used: 



log(ti(y,R))^log{y/\+4m- 

n foTy<R-l 
\ \(l-y + R) foTR-l<y<R+l 



(5) 







foiy > R + I 



The redshift interval from the observer to the source is divided 
into lens planes which are randomly populated by microlenses 
with masses and velocities following some arbitrary distribu- 
tion. All lens planes are evenly distributed in redshift between 
the source and observer, with a linear separation given by (dis- 
tance by light travel time): 



D' 



Ac 



Ho j 



dz 



il+z)^^Q(z) 



where 



Q{z) = Om(1 + zf - (Om +^a- 1)(1 + zf + 



(6) 



(7) 



Due to circular symmetry, it is sufficient to consider rectangular 
lens planes centered on the source with lenses moving in one 
dimension only. If we assume the lenses to move along the x- 
axis, there are two relevant boundaries Z, of the lens plane in 
the perpendicular direction, depending on whether the source 
is small or large. If Re represents the Einstein radius of the 
lenses, 



ly 



/ 27?Eymax for 7? < 1.2 

\ 2(RQsoiDd/D,) + RE) foiR > 1.2 



(8) 



where Z, in the case of R > 1 .2 corresponds to twice the size of 
the boundary of the outer non-zero magnification zone in (^. 
In the case of R < 1.2, S93 shows that y^ax = 15 is sufficient 
to ensure that the contribution to the magnification from lenses 
outside the considered region is negligible. 

In the case of an extended distribution of lens masses, ly 
should be set equal to whichever of the two expressions in (^ 
that produces the highest value when evaluated for the maxi- 
mum Re of any mass considered. 

The extent of the lens plane along the x-axis is derived by 
adding ly to a term which ensures that no lens originally outside 
the lens plane will move sufficiently close to the source during 
the time span At of the simulation to give a non-negligible con- 
tribution to the magnification (see S93): 

l_y = ly + 2Af(v'ob.s + max |Vcompactl), (9) 

where Zd is the redshift of the lens plane and Vobs and Vcompact are 
the observer and lens velocities, respectively, perpendicular to 
the line-of-sight. In the case of a Gaussian distribution of lens 
velocities, max Ivcompactl may be taken to be the 3cr-deviation. 

For the lens population, a constant comoving density is as- 
sumed. In the case where all compact objects share the same 
mass, the number of lenses populating a lens plane becomes: 



Pcrit,o(l + Zdy^compactD'f^Jxly 



compact 



(10) 



compact 



At each time step during the simulation, the observer, source 
and lenses are relocated according to their individual velocities 
and the total magnification calculated using (|l]). The resulting 
time-dependent magnification ;Utot(0 represents the simulated 
light curve. S93 may be consulted for further details on this 
method. 
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2.2. General assumptions 

In this paper, a A-dominated cosmology will be assumed, in 
which Q.M - 0.3, Q.a = 0.7 and Ho - 65 km/s/Mpc. The mean 
density inside the simulated volume is assumed to be equal to 
the mean density of the universe, which implies a homogeneity 
parameter (Kayser et al. 1997 )?] - 1 .0. All compact objects are 
furthermore considered to be randomly distributed. 

When demonstrating how the constraints on ficompact de- 
pend on different parameters, we will assume all compact 
lenses to have the same mass Mcompact. Unless stated otherwise, 
the velocities of compact objects perpendicular to the line-of- 
sight are assumed to be normally distributed with 2D velocity 
dispersions (Tv,compact = 400 km/s. 

Since the quasars used in this investigation are all located 
in the direction a - 21^28", 5 = -45°, a universal observer 
velocity will be used. Adopting the velocity of the Sun relative 
to the cosmic microwave background derived by Lineweaver 



0.18 



et al. ( 1996 1, the velocity of the observer perpendicular to the 



line-of-sight becomes 307.9 km/s. 



2.3. Tests 



The implementation of the algorithm has been subjected to sev- 
eral tests. It passes the test concerning the mean number of 
lenses with impact parameter y < 1 as outlined in S93. When 
the magnification formulas (^& ^ are modified to acco rd wit h 
the definition of magnification peaks used in Alexander ( 1995 ), 
the characteristic time scale of microlensing events as a func- 
tion of lens redshift also agrees well with the analytical expres- 
sions in that paper Finally, the standard deviation in magni- 
tudes Am around the mean magnification of the simulated light 
curves, calculated as 



Am 



m{ty - m{t) 



(11) 



where the deviation in magnitudes m{ f) from the mean is given 
by 



m{f) = 2.5 log//tot(f) - 2.5 logyUtot(f), 



(12) 



almost perfectly matches the analytical expression derived for 
large sources in Surpi et al. (2002), as demonstrated in Fig. |l]. 

For the current investigation, 30 lens planes have been used. 
The adequacy of this number has been verified by increasing 
the number of lens planes by a factor of two without any sig- 
nificant impact on the results. 

2.4. Flux conservation 

Magnification, as used here, is defined relative to the smooth 
Friedmann-Lemaitre universe, which is not a consistent treat- 
ment in terms of flux conservation. Even though this inconsis- 
tency has no effect on the variability of simulated light curves, 
the properly normalized scale of magnification becomes im- 
portant when simulating the effects of amplification bias. To 
force Ji^t - 1 when taking the average along all lines of sight 




Fig. 1. Analytical Ajn (Surpi et al. 2002 ) compared to Am de- 
rived from simulated light curves. At each quasar redshift, 
Zqso, the Am derived from simulations has been based on 50 
generated light curves, each spanning 2000 years and sampled 
at intervals of 1 year. Lines represent the analytical predic- 
tions for a A-dominated universe (Qm = 0.3, Qa = 0.7) with 
^^compact = 0.3 (soHd) and Ocompact - 0.15 (dashed). In both 
cases, we have assumed Rqso - 3 ■ lO'^ m. Circles and trian- 
gles indicate the corresponding results from simulations. 

in which the corrected magnification fi^°"{t) is calculated from 
the uncorrected yu™™'^(f) using 



Here, /con is a correction factor given by 

/con- = l/^CP(0 



(13) 



(14) 



where the average //""°"(f) is taken over all light curves gen- 
erated with randomized lens positions but otherwise identical 
input parameters. 

3. How to infer upper limits on the cosmological 
density of compact objects from the long-term 
variability of quasars 

3.1. General metliod 

Even though the observed optical variability of quasars on time 
scales of a few years (e.g. H ook et al. 1994 Veron & Hawkins 
Cristiani et al. 



1995 



1996) could be a combined effect of in- 



to a source, we foUow a recipe developed by Canizares (1982) 



trinsic (e.g. accretion disk instabilities, supernova explosions) 
and extrinsic (e.g. microlensing) variations, combining the two 
proposed mechanisms can however only increase the probabil- 
ity of flux variations. By assuming that all variability is due to 
microlensing and comparing the predictions from microlensing 
scenarios to the observed variation probabilities of quasar light 
curves, upper limits on the cosmological density of compact 
object may therefore be imposed. This technique was first im- 
plemented in S93 to constrain compact dark matter populations 
in the mass range 10"^ Mq-1 Mq for an EdS universe, using 
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the observational sample of Hawkins & Veron ( 1993 hereafter 
HV93). 

The method for constraining Qcompact outlined in S93 is 
based on the amplitudes, 6m, of quasar light curves. Here, the 
amplitude is defined as the difference between the minimum 
and maximum yearly magnitudes observed in a quasar within 
the duration of the monitoring programme. Each yearly magni- 
tude is furthermore the average of roughly four intrayear mea- 
surements. For the case of variations induced purely by mi- 
crolensing, this reduces to: 



6m = max(2.5 log/Ztot(0) - min(2.5 log//tot(0). 



(15) 



where the yearly magnitudes 2.5 log /itotCO for each integer year 
t are formed from four evenly spaced intrayear magnification 
data points. 

When comparing observed and synthetic quasar samples, a 
simple statistical analysis is used. The cumulative probability 
P(6m) predicted by microlensing simulations of finding ampli- 
tudes higher than 6m is used to calculate the expected number 
of objects in an observed sample of size A^sampie with amplitudes 
higher than 6m: 

< N(6m) A?sample^(<J'«). (16) 

The Poisson probability, P(N{6m)), that the number of ob- 
jects in the synthetic sample with amplitudes higher than 6m 
is smaller than, or equal to the observed number N{6m) then 
becomes 

-<N(Sm)> < N(6m) >'• 



P{N{6m)) = e- 



n=() 



(17) 



(18) 



The minimum of this probability 
P^in = minW5m)), 

6m 

essentially corresponds to the probability that the missing frac- 
tion of high amplitudes in the observed sample is simply due to 
the finite number of observed objects. Here, it will be used as a 
measure of the probability that the observed sample is compat- 
ible with a particular microlensing scenario. 

In principle, the statistical properties of synthetic samples 
depend on a large number of parameters (for the background 
cosmology Dm, ^^a> Hq and rj; for the lens population Qcompact. 

compact ^nd parameters describing the mass spectrum of com- 
pact objects; for the quasar population 0-^ qso and the distribu- 
tion of zqso in addition to parameters describing the luminosity 
function and the Lqso-^qso relation; for the observer the veloc- 
ity perpendicular to the line-of-sight). If most of these can be 
constrained by measurements or reasonable assumptions, up- 
per limits on ^compact may be imposed. In order to facilitate the 
visualization of the constraints imposed by this method, and 
still enable a direct comparison to the results of S93, a threshold 
value of Pn,in = 10% has been adopted to delimit allowed and 
rejected regions of the microlensing parameter space. When 
Pmm falls below this value, the microlensing scenario respon- 
sible is considered to be inconsistent with the observations and 
can be ruled out. A limit of P^in = 10% corresponds to the case 
where the microlensing scenario has predicted the presence of 
2.3 quasars with amplitudes higher than some arbitrary value 
in a sample, yet none is observed. 



3.2. Transition to the A-dominated universe and tine 
use of the large-source amplification formula 

What impact does the transition from an EdS (Qm - 10, 
Q.A = 0.0) to a A-dominated universe (Qm - 0.3, Qa - 0.7) 
have on the constraints derived? Since the A-dominated cos- 
mology implies a higher angular size distance out to a par- 
ticular redshift, more lenses will contribute to the light curve 
for a particular value of Qcompact- In Zackrisson et al. (2002) 
we showed that the constraints on Qcompact were significantly 
strengthened due to this effect. 

Since the approximation (Eh used in S93 underestimates the 



magnification in the limit of large sources, Surpi et al. (2002) 
predicted that the implementation of (@) should improve con- 



straints of S93 type. In Zackrisson et al. ( |2002[ ), it was however 
shown that in the range Mcompact = 10"'' Mq-1 Mq this correc- 
tion only had a modest impact on the upper limits on Qcompact 
in a part of parameter space (Mcompact = 10"'* M©, Rqso > 10'^ 
m) for which no interesting constraints could be inferred. 

As becomes evident when comparing the constraints from 
Zackrisson et al. (2002) to those of S93, there exists a system- 
atic difference between the cumulative probabilities P{6m) de- 
rived when analyzing the samples of HV93, even when the sim- 
ulations are matched as closely as possible given the published 
information. It seems that the model of S93 systematically pre- 
dicts somewhat higher probabilities of large amplitudes. For 
this reason, the constraints derived for the A-dominated cos- 
mology are actually only slightly stronger than those originally 
presented in S93. The origin of this discrepancy is not known, 
but could be due to the redshift resolution or the value of Hq 
used. 



3.3. Uncertainties in Rqso 

The microlensing magnification of a distant source is highly 
sensitive to the source size assumed. Even though the maxi- 
mum magnification from a single lens decreases with increas- 
ing source size, more lenses may at the same time contribute 



to the microlensing effect. As shown in Surpi et al. (2002), the 



variability amplitude does however drop as /?qso is increased 
for a fixed lens mass within the parameter space explored. 

The typical value of the parameter /?qso is far from well- 
determined, but certain constraints do exist. If we assume the 
UV-optical continuum of quasars to originate in the accretion 
disc surrounding the central black hole, we may expect /?qso to 
be at least a few times larger than the gravitational radius of the 
central mass, which for an 10** objects amounts to 1.5 ■ lO" 
m. 

Tight upper limits on Rqs,o are also available for a few well- 
studied quasars for which measurements of duration or magni- 
fication of microlensing events have been made. For the gravi- 



tational lens system QSO 0957H-561, Refsdal et al. ( |2000| ) have 
been able to place the constraint Rq^o < 6 ■ lO'^* m at a signif- 
icance level of 10%. For the same system Pelt et al. ( 1998| ) 
have argued that Rq^o ~ 3 ■ lO'^ m (within a factor of 3). 
Shalyapin (2001) has similarly inferred Rqso < 1 • 10'^ m for 
QSO 2237H-0305, with a most likely value at Rqso ~ 6 • 10'^ 
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m. It is however not clear how typical these values are of the 
quasar population as a whole. 

The constraints on ficompact quoted from S93 are based on 
the assumption that Rqso - 10'^ m. In the following, we will 
however investigate to what extent the constraints on Qcompact 
are affected when Rqso is allowed to take on higher and lower 
values inside the range Rqso - lO'^-lO'"^ m, which is an in- 
terval commonly considered in studies of quasar microlensing 
(e.g. Tadros et al. |1998[). 



and a = -3.41, yS = -1.58, M;(0) = 
-0.27, O* = 0.36 ■ 10-'' Mpc-^mag-i. 



-22.65, ki = 1.36, k2 



3.4. Building syntlietic samples 

In addition to the parameters that determine the variability of an 
individual quasar, distributions of vqso, Zqso, ^qso and Lqso 
also need to be considered when building the synthetic quasar 
samples used for comparison with the observed amplitude dis- 
tributions. 



3.4.1 . The velocity distribution of quasars 

Even though the clustering properties of quasars as a function 
of redshift are poorly determined, the tendency to cluster does 
not appear sufficiently strong to jeopardize a comparison with 
the kinematics of galaxies. We will therefore use galaxy data to 
estimate the velocity dispersion ct^ qso of quasars perpendicu- 
lar to the line-of-sight. 

In an investigation of synthetic galaxy catalogues from the 
Virgo consortium (Coil et al. 2001) one finds a typical line- 
of-sight velocity dispersion of galaxies within a 1 Mpc box at 
z = 1 of (Ty ~ 180+20 km/s. If large scale motions are included, 
the value increases to roughly 300 km/s. Empirical results from 
the Las Campanas survey (Baker et al. 2000), corresponding to 
low z, result in slightly lower values. Based on these numbers, 
and taking into account the fact that quasars are predominantly 
located at high z, we will adopt a universal quasar velocity dis- 
persion of cTy Qso = 300 km/s. 



3.4.2. Amplification bias 



A quasar that lacks sufficient intrinsic brightness to be included 
in a flux-limited sample may still be magnified through grav- 
itational microlensing to reach above the threshold for detec- 
tion. A flux-limited sample will therefore contain an enhanced 
fraction of highly magnified objects. The correlation between 
minimum magnification and amplitude pointed out by S93 im- 
plies that a flux-limited sample should also display an enhanced 
probability for large variations. When deriving P{6m) for a syn- 
thetic sample, the eff'ects of amplification bias should ffierefore 
be taken into account. 

In order to simulate the impact of amplification bias on the 
flux-limited samples of HV93 and H2000, we have assumed 
the optica l quasar luminosity function (LF) derived by Boyle et 
al. (2000) for a Q.u - 0.3, Qa - 0.7 cosmology: 

4)(Mb,z) - jq0.4[(q'+1)(Mb-M*(z))] Kf.m+DiMti-Mlim ^^'^^ 

where the redshift evolution is given by 

= ^b(O) - 2.5(^iz + k2Z^) (20) 



Using this LF, the probability for a quasar to have a certain 
Mb as a function of z may be derived. 

When building a synthetic sample subject to amplification 
bias, a large number of light curves is first generated with 
uniform zqso distribution within the redshift span of the ob- 
servational sample. For each of these, an absolute magnitude 
following the derived Mb probability is then randomly gener- 
ated. These Mb are converted into intrinsic apparent magni- 
tudes ;7iB,o assuming a power-law continuum with slope a - 0.5 
for all quasars when calculating the k-correction. The eff'ect of 
emission lines and deviations from a simple power-law con- 
tinuum on the k-correction has been neglected, since at least 
for 0.4 < zqso < 2.2 objects the a - 0.5 k-correction closely 
resembles the average k-correction derived from observations 
(Wisotzki|2000|). 



Since the observed quasar magnitudes are defined differ- 
ently in the HV93 and H2000 samples, the mB,o associated 
with each synthetic light curve is corrected for the effects of 
microlensing using different formulas depending on the obser- 
vational sample used for comparison; 



min(2.51og(//tot(0)) forHV93 
2.5 log(/ztot(0)) for H2000 



(21) 



where 2.5 log(//tot(0) are the yearly magnitudes defined in sec- 
tion 3.1. All objects fainter than the threshold used to define 
the observational sample are then rejected for efficiency rea- 
sons. This procedure is repeated until a sufficiently large num- 
ber of synthetic light curves have passed the threshold. Finally, 
the actual niB and zqso distributions of the observed samples 
are reproduced by dividing both synthetic and observational 
samples into identical bins and for each bin randomly select- 
ing a number of objects from the synthetic sample proportional 
to the relative content of the corresponding observational bin. 

This procedure of generating synthetic samples with simu- 
lated amplification bias differs from the one adopted by S93 in 
the assumed LF, the enforced flux conservation and the more 
realistic and zqso distributions of the final samples. A com- 
parison between f derived from synthetic samples with and 
without enforced observational flux distributions shows that the 
amplification bias produced by this method turns out to be quite 
modest, and most closely resembles the s - (weakest bias 
considered) scenario explored in S93. 

In simulating the amplification bias, we have assumed the 
observed LF to be identical to the intrinsic one, thereby ne- 
glecting the effects that lensing could have on the observed LF 
itself. Since ffie effect of lensing by compact objects as well as 
isothermal galaxies and clusters is to decrease the proportion of 
low-luminosity quasars compared to the intrinsic LF (e.g. Pei 



1995), we are in fact underestimating the amplification bias by 



adopting the observed LF. This ensures ffiat ffie upper limits on 
incompact derived are conservative. 
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3.4.3. The Lqso-^qso Relation 

In relating the intrinsic luminosity of a quasar to the size of the 
UV-optical continuum-emitting region, two different models 
have been considered. 

First, we have considered one in which Rqso is assumed 
unrelated to Lqso and treated as a free parameter in the range 
10'2-10''*m. 

We have also experimented with the thin accretion disk re- 
lation of Czerny et al. (1994) which, assuming a power-law 
continuum with index a - 0.5, may be written (SI units): 



Rqso = 3.37 ■ 10 



20 eff ieff 
c2(l 



(22) 



In this study, the effective wavelength Aes has been taken to 
be 4400 A since we are dealing with measurements in the B- 
band. Assuming Johnson filters and a - 0.5, L^^^ (J/s/m) is 
approximately related to Mb through 



= 6.77 ■ 10'* ■ lO^-/-^ \ 



(23) 



When generating the synthetic samples in the case where Lqso 
is assumed related to Rqso, ( ^2| ) is used in conjunction with 
( 23 ) to assign each quasar with a certain intrinsic Mb to one of 
the light curves with the closest discrete value of Rqso from the 
grid described in section 3.4.4. 



3.4.4. Parameter values and size of synthetic samples 

When generating the grid of synthetic samples used in the com- 
parison with the observations, all combinations of the discrete 
parameter values listed in Table 1 have been used, giving a to- 
tal of 150 simulation configurations. To ensure that the mag- 
nification probabilities derived from each sample is statisti- 
cally significant, a large number of light curves must be gener- 
ated for each such parameter combination. In this study, 35000 
light curves with a uniform redshift distribution in the range 
Zqso = 0.13-3.6 (Azqso = 0.01 between each bin) have been 
considered adequate. Simple tests show that this number results 
in Pmin derived from the final samples which are sufficiently 
stable for the present purpose, i.e. to establish the approximate 
border between rejected and allowed regions of the microlens- 
ing parameter space. 



Table 1. The set of discrete parameter values which define the 
grid of synthetic light curve samples used in this study. 



^compact • 

logMcompac, (Ma): 
log^Qso (m): 



0.05 

-4 

12 



0.1 

-3 

12.5 



0.15 

-2 

13 



0.2 
-1 

13.5 



0.25 



14 



0.30 



4. Upper limits on the cosmological density of 
compact objects 

4.1. Upper Limits on ^compact inferr ed fro m tine 
samples of Hawkins & Veron (|)995[) 



In HV93, two observational quasar samples are described: a 
FOCAP sample with a limiting magnitude of -21 and 
a bright sample with threshold otb = 18.5. All objects in 
these samples have been monitored for a period of 10 years. 
Following S93, these two samples have been combined into 
one, consisting of 117 objects in the redshift range zqso - 
0.29-3.23. In the context of deriving upper limits on ^compact 
from light curve amplitudes, this HV93 sample is however not 
ideal, since the variability is not explicitly expressed as an am- 
plitude. Instead, a related parameter s is used, which may only 
approximately be transformed into an amplitude through the re- 
lation 6m X s/lQ. Because of this complication, the constraints 
on Ocompact derived from the HV93 sample should be regarded 
as less certain than those inferred from the samples of H2000, 
and are only included here to facilitate a comparison of this 
study to that of S93. 

When calculating the synthetic 6m distribution, the time 
sampling must also be carefully considered. Even though the 
measurements of HV93 span 10 years, the s parameter is only 
derived from the 7 years when more than one plate was taken. 
This implies that the amplitudes of the synthetic light curves 
should be derived from 7 data points with the same separation 
as in HV93 (Hawkins 2002, private communication), not from 
the 1 1 evenly spaced data points assumed in S93. The effect of 
using too many yearly data points is to increase the probability 
of observing large variations, thereby making the constraints 
inferred on Qcompact too strong. When calculating the ampli- 
tudes, S93 furthermore neglects the averaging over intrayear 
magnitude variations described in section 3.1, thereby infer- 
ring to strong constraints on scenarios in which the typical time 
scale of variations is smaller than one year 

In Fig. ||we display the limits on fJcompact inferred from the 
HV93 sample when the more realistic time sampling described 
above is used and assuming that Rqso may be considered a 
free parameter As seen, the efficiency of this method to im- 
pose upper limits on ficompact is very sensitive to the value of 
Rqso assumed. For Rqso ^ 1 ■ 10'^ m, strong constraints of 
incompact ^ 0.05-0.1 may be imposed on compact objects in the 
mass interval Mcompact - 10"-' Mq-IO"^ Mq. At Rqso = 3-10'^ 
m however, the upper limits are significantly weaker and at 
^Qso - lO'** m no meaningful constraints may be imposed. 

In Fig. |]we also indicate the constraints on Ocompact when 



the Lqso-^qso relation ( |22| ) is used. The constraints are iden- 
tical to those produced for Rqso = 10'^ m, which is explained 
by the strong peak located at this value in the predicted Rqso- 
distribution. Fig. ^illustrates the 7?qso -distribution of the syn- 
thetic HV93 sample when Qcompact = 0.3 and Mcompact - 
10"^ Mq. The distributions for other parameter combinations 
only differ slightly due to the effect of amplification bias. 
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0.3 r 



0.25- 



0.2 



^ 0.15- 



0.1 - 



0.05- 




log M . (solar masses) 

compact 



Fig. 2. Upper limits on ^compact for five different lens masses 
inferred from a comparison with the combined samples of 
HV93 when Rqso is assumed unrelated to Lqso- The differ- 
ent lines represent Rqso = 10'"* m (thick solid), 3 ■ 10'"* m 
(thick dashed), 10'^ m (solid), 3 ■ 10'^ m (dashed) and 10'^ 
m (dash-dotted). Slight offsets from the parameter values of 
Table 1, for which the formal upper limits are derived, have 
been introduced to prevent overlapping lines. The maximum 
upper limit of Ocompact = 0.3 is here set by the background 
cosmology (Qm = 0.3, Qa - 0.7). 



4.2. Upper limits onQcompact inferred from tlie samples 
of Hawkins (|2000D 



We now shift our attention to the larger samples of H2000, 
which contain a total of 386 objects in the range zqso = 0.136- 
3.58, monitored for a period of 20 years. H2000 contains three 
subsamples: UVX, VAR and AMP. 

The UVX sample, which contains 184 objects in the red- 
shift range z - 0.242-2.21 is formed from a criterion of ul- 
traviolet excess, U - B < -0.2 and a limiting magnitude of 

— 21.5. 

The first variability-selected sample VAR, containing 298 
objects in the redshift range z - 0.136-3.58, has a limiting 
magnitude of -21 and an amplitude threshold of 6m > 
0.35. The second variability-selected sample AMP, containing 
66 objects in the redshift range z = 0.14-2.07, has a limiting 
magnitude of -21 and an amplitude threshold of 6m > 1.1. 

The constraints on Qcompact that may be derived from a com- 
parison of microlensing simulations to the statistical properties 
of these samples depend strongly on their selection criteria and 
redshift distributions. In order to find the most powerful upper 
limits, several combinations of the different subsamples have 
been tested. 

In the following, the light curves of the synthetic samples 
have been sampled at 21 evenly spaced yearly data points, each 
constituting an average of four intrayear magnification points, 
in fair agreement with the observational procedure (Hawkins 
2002, private communication). 

In Fig. ^ we present the constraints on Qcompact inferred 
from the UVX data set when assuming Rqso to be unrelated 



0.3 



0.25 



r, 0.2 



^ 0.15 



0.1 



0.05 




- - 3 ■ 10 m 



-1 



log M 



(solar masses) 



Fig. 3. Same as Fig. | for the UVX sample of H2000. 



to Lqso- The constraints are weaker than those produced from 
a comparison to HV93 (Fig. |]). The most important reason 
for this is the different observational selection criterion used, 
which prevents high redshift objects from entering the UVX 
sample. At z > 2.2 the Lyman forest enters the U-band, ren- 
dering the quasars red and out of reach of the U - B < -0.2 
criterion. 

It is reasonable to expect that a comparison to the 
variability-selected samples should result in more powerful up- 
per limits on Qcompact than those inferred from UVX, since 
these samples may extend to much higher redshifts, and the 
probabilities of high amplitudes increase with redshift in the 
microlensing paradigm, but not - as shown in H2000 - in the 
observations. The main obstacle of this approach lies in the 
amplitude thresholds, which overestimates the true variation 
probability by rejecting low-amplitude objects. If the amplitude 
threshold is too high, no meaningful constraints on Qcompact can 
be imposed. 

One may also combine UVX and variability selected sam- 
ples to form a data set which both extends to high redshifts 
and is less Ukely to substantially overestimate the variation 
probability of quasars. In Fig. Q we indicate the constraints on 
Qcompact inferred from the combined UVX-hVAR sample when 
Rqso is assumed unrelated to Lqso- Due to the extension to 
higher redshifts, the upper limits are significantly stronger than 
those inferred from UVX. The upper limits derived from the 
VAR sample alone are somewhat weaker than those presented 
in Fig. |. 

The AMP sample is not useful for these purposes, both be- 
cause of its very high amplitude threshold and its limited red- 
shift span. 

Since the probabilities of high amplitudes increase with 
redshift in the microlensing scheme, a possible way to improve 
the efficiency of this method even further may however be to 
conduct a comparison between a high-redshift subset of the 
VAR sample and the microlensing simulations. 

H2000 shows that the higher redshift objects in the VAR 
sample display less variability than their low redshift coun- 
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terparts. This is attributed to the fact that the higher redshift 
quasars also have higher luminosities and the hypothesis that 
Rqso increases with Lqso- However, by stretching the assump- 
tion that these two parameters are unrelated, the fact that the 
high redshift quasars of VAR show smaller amplitudes than 
their low redshift counterparts should not hamper the exclu- 
sive use of high redshift objects in the comparison. Fig. |] 
shows the substantially improved constraints when all objects 
with z < 1 .5 have been removed from the VAR sample. This 
leaves 144 objects, which still provides better statistics than the 
samples of HV93. The constraints inferred this way are very 
strong: Qcompact < 0.05 for Mcompact = 10"^ Mq - 10"' Mq and 

f^compact < 0.1 for Mcompact = 1 Mq aS long aS /?QSO < lO'^ 

m. Even at Rqso = 3 ■ lO'-' m competitive constraints may be 
derived for Mcompact = 10"^ Mq - 10"' Mq. 

Even though it may be possible to construct a Lqso-^qso 
relation which mimics the inverse coiTelation between lumi- 
nosity and amplitude seen in the observed samples, while al- 
lowing essentially no amplitude dependence on redshift for 
z > 0.5, the Lqso-^qso relation explored here does not have 
this property. In Fig. ^ we indicate the constraints on Qcompact 
derived from a comparison to the z > 1 .5 part of VAR in the 
case when the Lqso-^qso relation ( p2| ) is used. In Fig. ^ we 
show the corresponding 7?Qso-distribution of the sample gener- 
ated when ^compact - 0.3 and Mcompact - 10"** Mq. Despite the 
higher redshift, the typical Rqso value actually lies somewhat 
lower than that predicted for the HV93 sample. The 1ow-7?qso 
tail not seen in the HV93 distribution is a result of the stronger 
amplification bias induced at high redshift, where the average 
intrinsic quasar luminosity lies further below the limiting lumi- 
nosity of the sample. 



5. Discussion 

As indicated in Fig. || the method of using the long-term vari- 
ability of quasars provides excellent possibilities to constrain 
Qcompact, especially with the use of high-redshift quasar sam- 
ples. For certain combinations of Mcompact and Rqso, Pmin falls 



far below the threshold value of 10% even at the lowest values 
of Qcompact Considered, indicating that the method could possi- 
bly be refined to constrain cosmological densities even below 
incompact - 0.05. Similarly, we note that it is not impossible that 
constraints on Mcompact = 10"^ Mq may also be imposed. Due 
to the very CPU-demanding simulations necessary at such low 
masses (many lenses contributing to the light curve), we have 
however not investigated this possibility here. 

Despite the improvement of the upper limits imposed here 
compared to those derived in S93, several uncertainties still 
prevent us from drawing definite conclusions about Qcompact in 
the mass range explored. 

As already shown, the constraints are highly sensitive to 
the typical value of the Rqso parameter. Even though the few 
available measurements of Rqso indicate values inside the re- 
gion for which powerful constraints may be derived, extrapo- 
lation of these results to the whole quasar population could be 
hazardous. 

Other parameter uncertainties also deserve consideration. 
The upper limits on Qcompact presented so far rely on the as- 
sumption that (Tv.compact - 400 km/s. This parameter essentially 
describes the preferred whereabouts of the compact objects. If 
the compact objects cluster only weakly, i.e. are mainly located 
in the field, a value of o-y compact = 200 km/s should be more 
appropriate. If they cluster more strongly, i.e. are mainly dis- 
tributed inside rich galaxy clusters, much higher values could 
in principle be considered. An estimate of the upper limit to 
this value may be obtained from model predictions by Sheth et 
al. ( 2001 ) under the assumption that the lenses dynamically be- 
have similar to cold dark matter particles. Based on the values 
they derive for different spatial scales we choose as an upper 
limit (Tv^corapact - 600 km/s. 

In the context of microlensing, the effect of varying 
o"v,compact IS to alter the characteristic time scale of source cross- 
ing and the width of the corresponding magnification peak. 
Higher velocities imply shorter time scales and, in principle, an 
increased probability of detecting high amplitudes since more 
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Fig. 6. Upper limits on Qcompact inferred for different Mcompact 
from the HV93 sample (thick solid) and the z > 1.5 part of 



VAR (thick dashed) when the Czerny et al. ( 1994 ) Lqso-^qso 
relation is applied. 



compact object may pass the source during the time span of the 
observational programme. The impact of this effect is however 
weakened by the finite sampling rate of the light curves. 

By expanding the grid of synthetic samples defined by 
Table 1 to other (Tv^compact at Rqso - 10'"' m, we have inves- 
tigated in what way the upper limits on Qcompact inferred from 
the z > 1 .5 subset of the VAR sample are affected by crv,compact 
variations in the range 0-^ compact - 200-600 km/s. The effect 
turns out to be very undramatic and only relaxes our upper lim- 



its from H, 



compact 



0.1 to 0.15 at M, 



compact 



1 Mo in the case of 



o"v,compact - 200 km/s. All other constraints at this source size 
are unaffected. 

Additional uncertainties stem from approximations inher- 
ent in the microlensing model used. The model assumes no 
shear and a circular source with uniform brightness. The inclu- 
sion of shear terms is expected to lower Am (Refsdal & Stabell 



1997), thereby making the constraints on Qcompact weaker. Due 



to the lack of analytical approximations for fj.(y,R) in the case 
of both non-zero shear and more realistic (e.g. Gaussian) source 
brightness profiles, the quantitative impact of such features can 
however presently only be tested with more time-consuming al- 
gorithms like backwards ray-tracing. The model also assumes 
all compact objects to be randomly distributed, i.e. that cor- 
relations between lenses are unimportant. However, for the 
most extreme combination of lens parameters considered in 
this study, the simulations include compact objects initially as 
far away as two light years from each other in the direction 
perpendicular to the line-of-sight. This means that compact ob- 
jects which are bound together on scales smaller than this (e.g. 
in clusters or binary systems) will have coiTelated positions in 
the lens plane. Simple tests show that correlations of this kind 
generally increase the probability of high amplitudes, which 
indicates that the upper limits derived from random lens distri- 
butions should be conservative. 
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Fig. 7. Examples of the Rqs.o distributions of the synthetic 
samples generated to match the HV93 sample (solid) and the 
z > 1.5 subset of VAR (dashed), when the Czerny et al. 



( 1994 ) Lqso-^qso relation is applied. In both cases we as- 
sume Qcompact = 0.3, Mcompact = lO"'* Mq and usc identical 
bin widths when plotting the distributions. 

6. Conclusions 

In this paper, the analysis of the upper limits on the cosmo- 
logical density of dark matter in the form of compact objects 
inferred from the long-term variability of quasars has been im- 
proved compared to that of S93 in several ways: 

- By generaUzation to the Qm - 0.3, Qa = 0.7 universe; 

- By implementation of a magnification formula more suit- 
able for microlensing of large sources; 

- By implementation of an improved method of generating 
synthetic quasar samples which matches the zqso and niQ 
distributions of the observational data sets; 

- By considering uncertainties in the typical value of Rqso 
and cr v,compact5 

- By taking the intrayear magnitude averaging of the ob- 
servational procedure into account when deriving the light 
curve amplitudes; 

- By comparing results from numerical simulations to larger 
observational samples with longer time span of quasar 
monitoring. 

We show that by using a high redshift subset of variability- 
selected samples, previous limits may be substantially im- 
proved. We infer Qcompact < 0.05 for Mcompact = 10"^ Mg- 

10"' M0 as long as typically Rqso < 1 ■ 10'^ m at z > 1 .5. The 
few existing measurements of Rqso are however not sufficient 
to ensure that this condition is fulfilled. At Rqso - 3 ■ 10'^ 
m the Qcompact ^ 0.05 constraint only holds for lenses in the 
mass range Mcompact - 10"^ Mq-IO"' Mq, and at Rqso - 10'"' 
m the only meaningful constraints are Qcompact 5 0.15 for 

Mcompact = 10"' Mo-1 Mq. 

When the strong /?qso -dependence of the upper limits 
is compounded with additional uncertainties discussed under 
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Section 5, we are forced to conclude that the method of us- 
ing the long-term variability of quasars to place upper limits 
on Ocompact cauuot be used to reliably rule out any cosmolog- 
ically significant populations of compact objects in the stellar 
to planetary mass range at the present time. This situation may 
of course change as more and better measurements of /?qso be- 
come available. 

Acknowledgements. We greatfuUy thank Patrik GuUerstrom for his 
contributions to the programming of the microlensing code. 

References 

Alam, J., Raha, S., & Sinha, B. 1999, ApJ 513, 572 
Alexander, T. 1995, MNRAS 274, 909 
Baker, J.E., & Davies, M. 2000, ApJ 536, 1 12 
Boyle, B.J., Shanks, T, Croom, S.M., et al. 2000, MNRAS 317, 1014 
Buries, S., Nollett, K., & Turner, M. 2001, ApJ 552, LI 
Canizares, C.R. 1982, ApJ 263, 508 
Carr, B.J., & Sakellariadou, M. 1999, ApJ 516, 195 
Coil, A.L., Davis, M., & Szapudi, L 2001, PASP 113, 1312 
Cristiani, S., Trentini, S., La Franca, R, et al. 1996, A&A 306, 395 
Czemy, B., Jaroszynski, M., & Czerny, M. 1994, MNRAS 268, 135 
Dalcanton, J., Canizares, C, Granados, A., & Steidel, C. 1994, ApJ 
424, 550 

Hawking, S. 1971, MNRAS 152, 75 
Hawkins, M.R.S. 1996, MNRAS 278, 787 
Hawkins, M.R.S. 2000, A&AS 143, 465 (H2000) 
Hawkins, M.R.S., & Veron 1993, MNRAS 260, 202 (HV93) 
Hook, I.M., McHahon, R.G., Boyle, B.J., Irwin,& M.J. 1994, MNRAS 
268, 305 

Kayser, R., Helbig, R, & Schramm, T. 1997, A&A 318, 680 

Lineweaver, C.H., Tenorio, L., Smoot, G.F., et al. 1996, ApJ 470, 38L 

Membrado, M. 1998, MNRAS 296, 21 

Mohapatra, R.N. 1999, Physics Letters B 462, 302 

Ostriker, J., & Vietri, M. 1983, ApJ 267, 488 

Pei, Y.C. 1993, ApJ 403,7 

Pei, Y.C. 1995, ApJ 440, 485 

Peh, J., Schild, R., Refsdal, S., & Stabell, R. 1998, A&A 336, 829 
Refsdal, S., & Stabell, R. 1997, A&A 325, 877 
Refsdal, S., Stabell, R., Pelt, K., & Schild, R. 2000, A&A 360, 10 
Rauch, K.P, & Blandford, R.D. 1991, ApJ 381, L39 
Sheth, R. K., Diaferio, A., Hui, L., & Scoccimarro, R. 2001, MNRAS 
326, 463 

Schneider, P 1993, A&A 279, 1 (S93) 
Shalyapin, V.N. 2001, Astronomy Letters 27, 150 
Surpi, G., Refsdal, S., & Helbig, P 2002, submitted to A&A 
Tadros, H., Warren, S., & Hewett, P. 1998 New Astronomy Reviews 
42, 115 

Veron, P, & Hawkins, M.R.S. 1995, A&A 296, 665 
Wisotzki, L. 2000, A&A 353, 861 

Zackrisson, E., Bergvall, N., & Helbig, P. 2002, in Lighthouses of the 
Universe: The Most Luminous Celestial Objects and Their Use 
for Cosmology, ed. Gilfanov, M., Sunyaev, R., & Churazov, E. 
(Springer), 550 



